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Abstract 

Here, dark energy is obtained using dual roles of the Ricci scalar (as 
a physical field as well as geometry). Dark energy density, obtained 
in this model, mimics phantom and the derived Friedmann equation 
contains a term Pj e /2A with pdc being the dark energy density and A, 
called as cosmic tension. It is like brane-gravity inspired Friedmann 
equation, which arises here without using the brane-gravity theory. It 
is found that acceleration is a transient phenomenon for A < 0, but 
for A > accelerated expansion is found to encounter the big-rip prob- 
lem. It is shown that this problem can be avoided if the dark energy 
behaves as barotropic fluid and generalized Chaplygin gas simultane- 
ously. Moreover, time for transition (from deceleration to acceleration 
of the universe) is derived as a function of equation of state parameter 
Wde = Pde/Pdc with p de being the pressure for dark energy fluid. 

Keywords : Dual role of the Ricci scalar; phantom dark energy;transient 
acceleration; Chaplygin gas and avoidance of big-rip singularity. 
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1. Introduction 

Astronomical observations, during the last few years, provide com- 
pelling evidences in favor of accelerating universe at present, which is 
caused by dominance of dark energy (DE) [1, 2, 3, 4, 5, 6, 7, 8]. Obser- 
vation of 16 Type la supernovae (SNe la) by Hubble Space Telescope 
further modifies these results and shows evidence for cosmic deceler- 
ation preceding acceleration in the late universe [9]. So, DE has a 
significant role in cosmic dynamics of the current universe. The sim- 
plest candidate for DE is supposed to be the cosmological constant A 
which is very high in the early universe, but there is no mechanism to 
bring it to present value without fine-tuning. Alternatively, to explain 
its decay from a very high value, in the early universe, to its present 
extremely small value, many models [10, 11] were suggested in which 
A is envisaged as a slowly varying function of cosmic time. Apart from 
dynamical A as DE, other models are fluid-dynamics models (where 
barotropic fluid is its source with or without dissipative pressure) , 
Chaplygin gas(GC) and generalized Chaplygin gas (GCG)models [12, 
13, 14, 15, 16, 17, 18, 19]. In field-theoretic models, the most natural 
ones are models, where DE is caused by scalars. These are quintessence 
models [20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32], k-essence mod- 
els [33, 34], tachyon models [35, 36, 37, 38, 39, 40, 41, 42, 43, 44] and 
phantom models [19, 45, 46, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57]. In 
these models, non-gravitational lagrangian density for exotic matter 
giving DE is added to Einstein-Hilbert term in the action. Recently, 
in a different approach, non-gravitational term is replaced by gravita- 
tional term being non-linear in Ricci scalar R, which stems modified 
gravity [58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 
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75, 76, 77, 78, 79, 80, 81, 82, 83, 84]. Thus, in these theories, non-linear 
term of R is introduced as DE lagrangian. 

In all these models, a non-gravitational or gravitational term for 
DE is added in the theory a priori. In this sense, these models are 
phenomenological. Here also, non-linear term of R is added to Einstein- 
Hilbert term, but contrary to [58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 

69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84], in this 
model, non-linear term of R is not used as lagrangian for DE. Rather, 
DE density emerges spontaneously, which makes it different from earlier 
models of gravitational DE [58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 

70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84]. In the present 
model , DE emerges as combined effect of linear as well as non-linear 
term of R manifesting dual role of the Ricci scalar as a physical field as 
well as geometry. Thus, here, DE density is derived from the modified 
gravity. So, this approach of getting gravitational DE is different from 
the approach of [58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 
73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84], where energy-momentum 
tensor of DE is obtained from non-linear term of curvature taken as 
DE lagrangian. Moreover, in what follows, it is found that dark matter 
also emerges due to non-linear term of curvature scalar. 

Natural units (h — c — 1) are used here with GeV as the fundamental 
unit, where h and c have their usual meaning. In this unit, it is found 
that lGeV" 1 = 6.58 x lQ- 25 sec. 



2. Gravitational action, modified Friedmann equation and 
cosmic cconsequences 
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The action for higher-derivative gravity is taken as 

S„ 



where R is the Ricci scalar curvature, G = Mp 2 {M P = 10 19 GeV is the 
Planck mass). Moreover, a is a coupling constant having dimension 
(mass) 2r . In what follows, it is shown that r = 1 + w de , where w de = 
Pdc/Pdc < — 1/3 with Pde(Pde) being pressure (energy density) for DE 
fluid. So, (2 — r) > showing that instability problem does not arise 
[58, 59, 60]. 

The action (1) yields gravitational field equations 
^(i^-^i*)-a[(2-r){V M V„i^ 

-\g^- r) \ = o (2) 

using the condition 5S g /5g^ v = 0. Here, denotes covariant derivative 
and the operator □ is given as 

1 d ( , ,.„ d 



□ 



^gdx^V dx h 
with /i, v — 0, 1, 2, 3 and g^ u as metric tensor components. 

Taking trace of eqs. (2) and doing some manipulations, it is obtained 
that 

r i?( 1+r ) i 

rR 2 } (4) 



OR- 4v"i2V„i2 = 



16nGa 



R 3(2-r)(l-r) 
with a 7^ to avoid the ghost problem. 

Experimental evidences support spatially homogeneous flat model of 

the universe [85, 86, 87, 88, 89]. So, the line-element, giving geometry 

of the universe, is taken as 

dS 2 = dt 2 - a 2 (t) [dx 2 + dy 2 + dz 2 } (5) 

with a(t) as the scale factor. 
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In the homogeneous space-time, given by eq.(5), eq.(4) is obtained 



as 



R + ^R-^ 
a R 



1 



r R {1+r) 

167iGa 



3(2 -r)(l -r) " * ' 
In most of the situations, for example, radiation model, matter- 
dominated model, and accelerated models, we have a(t) as a power- 
law solution yielding R as the power-law function of a(t) . So, it is 
reasonable to take R as 

A 



R 



(7) 



with n > being a real number and A being a constant with mass 
dimension 2. 

R, given by eq.(7), satisfies eq.(6), if 



1 



- + (2 + rn - n) - = — 
a V ; Va/ 3(2-r)(l-r) 



a- nr rA 

1 a 

87iGna n 



(8) 



Eq.(8) integrates to 



a\ ■ 



C 



+ 



1 



a) a 2[3-n(i-r)] 3(2-r)(l-r)L 87rGria(-rar + 2 [3 - n(l - r)])A r 
+ ^ . J A aH (9) 



ra(-n + 2[3-n(l-r)]) J 
with C being an integration constant. Setting n = 3 eq.(9) yields the 

modified Friedmann equation (MFE) 



/a\ 2 87rG / 



with 



8ttG 



"Pde = 



2167rGar(2-r)(l - r) 



(10) 



;n) 



and 



8ttG 

— Pm = 



27(2-r)(l-r)(-l + 2r)a 3 



(12) 



6 
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X = - =-. (13) 

48nGC{72nGar( y 2 - r)(l - r)} 2 

p m , given by eq.(12), has the form of matter (pressureless ideal fluid) 
density emerging from the gravitational sector. So, it is termed as 
density for cold dark matter. Moreover pd e , given by eq.(ll), is also 
inspired by modified gravity, so it is called DE density. 

Conservation equation for DE 

Pde + 3-p de (l + w dc ) = 0, (14a) 
a 

where equation of state (EOS) parameter 

W dc =Pdc/pde (146) 

with pde and pd c being isotropic pressure and density for DE respec- 
tively. 

Eq.(14a) yields 

r = (l + w dc ). (15) 
Connecting eqs.(12) and (15), it is obtained that p m > 0, if 

(1 + W dc ) 



W de (l - W dc )(l + 2w de ) 

The inequality (16) is possible for 



> (16) 



(l + w dc )<0 (17a) 

as DE characterized by Wd c < —1/3. It means that DE, obtained here, 
mimics phantom. Thus, eq.(14) and (16) yield 

(l-r) = -w dc >0. (176) 

WMAP data give the present value of pd c to be p° de ) = 0.73p cr . and p m 
to be p° n = 0.23p cr ., where p cr . = 3Hq/8itG with H = lOOakm/secMpc = 
2.33 x l(T 42 GeV and h = 0.68 (having maximum likelihood) [90, 91]. 
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Using these values, constants a and A in eqs.(ll) and (12) are evaluated 

as 

2167rGa(l - w de )| - w dc (l + w de )| = A r H 2 af +Wdcl (18) 

and 

27(l-w de )|-w de (l + w dc )| °' K J 

where a = a(t ) with t = 13.7Gyr = 6.6 x 10 41 GeV _1 being the 
present age of the universe. Here onwards, ao is noemalized as 

a = 1. (20) 

Incorporating eqs. (15)-(20) in eqs.(10)-(13), it is obtained that 

( ^^ Wl + t} + Pm] (21o) 

with 

Pdc = 0.73p cr a 3|1+Wdcl , (216) 



p m = 0.23p cr a" 3 , (21c) 

and 

rl6 -i — i 
A= — ttGC (21d) 

- O 

It is interesting to see that MFE (10) contains a term p^ c /2A, which 
is non-linear in p dc . As a reference, it is right to mention that brane- 
gravity inspired Friedmann equation also has a non-linear term like 
p 2 /2Xb with p being the energy density and A& being the brane-tension 
[92]. Here, MFE (10) arises from the modified gravity without us- 
ing any prescription for the brane-gravity and A is termed as cosmic 
tension. 
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In eq.(10), energy density is a physical concept, which is caused by 
the scalar curvature R. It shows role of R as a physical field also. So, it 
is reasonable to remark that R has dual role here (i) as a physical field 
and (ii) as a geometrical field [93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 
103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116]. 

A, given by (21d), can be positive as well as negative. So (21a) is 
re-written as 

with (+) sign for A > and (— ) sign for A < 0. In case, A < 0, effect of 
phantom DE vanishes when pd e grows to 2|A| in future and accelerated 
expansion will come to an end. 

Now, p m > pdc when red- shift z is given as 

/73\ 1 /3|w dc | 

z > y - 1. (23a) 

Pm < Pdc for 

/73\i/3|w dc | 

z < ( — - 1. (236) 



,23, 

Inequalities (23a) and (23b) show that transition from p m > pd c to 
Pm < Pdc takes place at 

Observations yield —1.22 < w dc < —1 [9] for phantom DE. Using 
these values of w dc , (24) yields 

0.37 < z* < 0.46, (25) 

which is supported by the range 0.33 < z* < 0.59 given by 16 Type 
Supernova observations [9]. 
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When z > z*, p m > pd c and pd e /2|A| < 1, so (22) reduces to 

0' = 0.28^(5)', ,26) 

using (21c) for p m . 
(26) integrates to 



/ 1 \ 3 / 2 

a{t) = a d 1 + 0.72^o ( — ) (* - *d) 



2/3 



(27) 



with ad = a(td) being a constant. a(t), given by (27), shows deceleration. 
Further, for z < z^.pm < pd c , so (22) reduces to 

Pdc 



= 0.73# V |1+Wdcl 



1 ± 



2IAIJ' 



Case 1. When A < 

In this case, (28) looks like 



ay 
a) 



0.73# V |1+Wdcl 



1 - 



Pdc 

2IAI 



which integrates to 



(28) 



(29) 



a(t) = 



^ + { a;3|1+Wdcl/2 V 1 



0.73p cr a 



3|l+w de | 



2|A| 

2n -l/3|l+w de | 



(30) 



giving accelerated expansion, which will continue tiil the time t e such 
that 

P(Q = 2|A|. 

implying 

0.73p cr af +Wdcl = 2|A|. (31) 
Connecting (30) and (31), it is obtained that 



3 1 1 + w dc |if 
giving period of acceleration. 



;i + ^) 3 i i + w dd _ 



0.73p c 
2IAI 



(32) 
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The time t* for transition from deceleration to acceleration is ob- 
tained as 



t* — CO" 



3|l+Wde|iW07T3 



;i + z„) 3|1+Wde| 



0.73p c 



21 A| 



0.73p c 
2|A| 
(33) 



For t > t e , pdc > 2 1 A | = p(t e ) and (22) takes the form 



a\ 2 87rG 
a 



- Pdc 



Pdc 

21 A I 



-1>+P, 



(34) 



It is interesting to see that (34) yields a — 0, at some future time 
t = t m ,when we have 

^ „, d \ „, (f \ 

(35) 



Pm(^m) _ Pdc(^m) ^ _ Pdc(^m) _ ^ 



Pde(^m) 2 1 A| Pdc(^c) 

and a(t) grows to its maximum a m satisfying the condition 



„3(l+w do ) _ 3(l+w de ) , f§ 3(l+2w do ) 
u e ~~ "m ^ 73 



(36) 



for (1 + w dc ) < 0. 

During the period (t m — t c ), (34) is effectively like (26) yielding de- 
celerated expansion, driven by matter, as 



a(t) = a c 1 + ^rVM3(t - t e ) 



2a e 3/2 



2/3 



Case 2. When A > 

In this case, (28) looks like 



-) 2 = 0.73tf Vl 1+w -' 
a) 



Pdc 



which integrates to 



(37) 



(38) 



a(t) = 



0-73p c 
21 A! 



+ 



■3|l+W de |/2, 



0.73p cr af +Wdc 



1 + 



21 A| 
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3|1+W dc l -I 2-1 -l/3|l+w de 



2 

showing a — > oo, when 



-fZo>/0J3(*-**)} J (39) 



t -> t s = U + [1.095|1 + w dc |# o r 1 x 



. (1 + ^ + 2^ ± ^]. (40 ) 

Moreover, as t — > t s ,Pdc — ► oo and p dc — > —oo showing the big-rip 
problem. 

3. Avoidance of big-rip singularity 

In what follows, it is shown that the "big-rip singularity", obtained 
above for positive cosmic tension A can be avoided if DE (obtained here) 
as a barotropic fluid obeying (14b) as well as generalized Chaplygin gas 
(GCG) with EOS given as [12, 13, 14, 14, 16, 17, 18, 19] 

Pde = 0— , (41) 

Pde 

where < (3 < 1 for GCG , /3 = 1 for CG and M is a constant. 

Due to negative pressure and being only fluid with supersymmetric 
generalization , GC and GCG have been strong candidates for DE. 
But experimental results support GCG [13, 14, 14, 16, 17, 18]. So, 
GCG model is preferred here. Avoidance of "big-rip singularity" , using 
double behaviour of DE fluid as being done here, was explored in [19] 
for non-gravitational phantom DE. 

(14b) and (41) yield 

w d c = -( — ) • (42) 

VPde/ 

This equation shows dependence of w dc on p de varying with a(t). Now, 
using eq.(42) for p dc with variable w de , conservation equation (14a) 



12 

looks like 
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W dc = 3-W dc (l + W de ), 



which integrates to 



w dc 



w. 





de_„3(l+/3) 



1 + w dc 1 + w[} c 

with w£ e = w dc (a ). 

Connecting (42) and (44) 



(43) 



(44) 



Pdc = M 



1 _ ( 1 +<) a -3(l+/3) 



W. 





de 



l/(l+/3) 



where 



M = 0.73p cr |w^ c 



|l/(l+/3) 



being obtained from (42). 

Now, for A > and p m < p dc , (22) takes the form 



a\2 8ttG 



"Pdc 



1 + 



Pdc 

21 A| 



(45) 



(46) 



(47) 



Connecting (45) and (47) and making approximations, it is obtained 
that 



a 8ttG 



a y 3 " V 2 
which integrates to 



1 - 



(l+Wde) n -3(l+/3) 



(48) 



a(f) 



+ 



{(1 + ^)" 



3(l+/3) _ 



(j + Zde) 



X 



3(1+/?) 



1/3(1+0) 



which is valid for 



(l + ^- 3(1+ ^>^-^ 



2w^ e 



(49) 



(50) 



(49) is free from any finite time singularity listed in [117]. Using (49) 
in (45) and (42), it is obtained that p dc — > M and w dc — > — 1. 
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4. Summary 

In summary, it is interesting to see that phantom like dark energy 
emerges spontaneously from the gravitational sector and the Fried- 
mann equation(giving cosmic dyanamics) contains a term p^ e /2A with 
Pde being the dark energy density and A being the cosmic tension . It is 
found that time of transition from deceleration to acceleration depends 
on Wde- Here, two cases arise (i) when cosmic tension A < and (ii) 
when cosmic tension A > 0. In case of negative cosmic tension, late 
cosmic acceleration begins from the transition time t* and stops at a 
finite time t s in future. It shows that the negative cosmic tension puts 
a brake on accelerated expansion at t — t s . Later on, universe decel- 
erates driven by matter showing another transition from acceleration 
to deceleration. Decelerated expansion contiues till a(t) acquires its 
maximum at t — t m > t s and the ratio p m (t m ) / pdc(t m ) is given by (35). 

When cosmic tension is positive, universe accelerates from the time 
of transition and ends up into 'big-rip'singularity in future. It is shown 
that this problem can be avoided if dark energy behaves as a barotropic 
fluid and GCG simultaneously. This type of approach was adapted in 
the case of non-gravitational phantom led cosmological model in [19] 
using classical mechanics. It is in contrast to approach in [118, 119, 
120], where it is suggested that "sudden future singularity" may be 
avoided using quantum gravity effects. 
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